Asymptotic Theory for the Probability Density Functions in Burgers Turbulence 
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A rigorous study is carried out for the randomly forced Burgers equation in the inviscid limit. No 
closure approximations are made. Instead the probability density functions of velocity and velocity 
gradient are related to the statistics of quantities defined along the shocks. This method allows one 
to compute the anomalies, as well as asymptotics for the structure functions and the probability 
density functions. It is shown that the left tail for the probability density function of the velocity 
gradient has to decay faster than |£| -3 . A further argument confirms the prediction of E et al. 
[Phys. Rev. Lett. 78, 1904 (1997)] that it should decay as |£| 



-7/2 



PACS numbers: 47.27.Gs, 05.40.-a, 02.50.Ey 

In this Letter, we focus on statistical properties of so- 
lutions of the randomly forced Burgers equation 



(1) 



where / is a zero-mean, statistically homogeneous, whitc- 
in-time Gaussian process with covariance 



(f(x,t)f(y,s))=2B(x-y)6(t-s), 



(2) 



where B(x) is smooth. We are particularly interested 
in the probability density function (pdf) of the velocity 
gradient £(x, t) = u x (x,t), since it depends heavily on 
the intermittent events created by the shocks. Assuming 
statistical homogeneity, and letting Q(£; t) be the pdf of 
£(x,i), it can be shown that Q satisfies 



Qt = CQ + (C 2 Q) f + BiQ« - v((&w\S)Q) 



(3) 



where B\ = —B xx (0). (£ X x\Q is the ensemble-average of 
£, xx conditional on £. The explicit form of this term is 
unknown, leaving (^) unclosed. There have been several 
proposals on how to approximately evaluate the quantity 



F{^t) = -\\mv{^ xx \i)Q) 



(4) 



At steady state, they all lead to an asymptotic expression 
of the form 



f c-\i\- a 

Q ~ 1 C + fe-Z 3 /^ 



as £ — > — oo, 
as £ — > +oo, 



(5) 



for Q, but with a variety of values for the exponents a and 
(3 (here the C±'s are numerical constants). By invoking 
the operator product expansion, Polyakov jjj suggested 
that F = aQ + b£Q, with a = and b = -1/2. This 
leads to a = 5/2 and (3 = 1/2. Boldyrev Q considered 
the same closure with — 1 < b < 0, which gives 2 < a < 3 
and (3 = 1 + 6. The instanton analysis ||,|| predicts 
the right tail of Q without giving a precise value for (3, 
but has not given any specific prediction for the left tail. 
E et al. H made a geometrical evaluation of the effect of 
F, based on the observation that large negative gradients 



are generated near shock creation. Their analysis gives 
a rigorous upper-bound for a: a < 7/2. In it was 
claimed that this bound is actually reached, i.e., a = 7/2. 
Finally Gotoh and Kraichnan j(| argued that the viscous 
term is negligible to leading order for large |£|, i.e. F « 

1 /3 

for £| ^> B^ . This approximation leads to a = 3 and 
[3 = 1. For other approaches, see e.g. fj|p|. In this letter 
we proceed at an exact evaluation of Mj and we prove 
that a has to be strictly larger than 3 (a result which 
does not require that steady state be reached) . At steady 
state, we prove that (3=1 and we give an argument 
which supports strongly the prediction of ||, namely, 
a = 7/2. 

To begin with, let us remark that it is established in 
the mathematics literature that the inviscid limit 



u°(x,t) 



lim u(x, t), 



(6) 



exists for almost all (x,t). Since u° will in general develop 
shocks, say, at x = y, we may have m° cx S(x — y), and 
one cannot simply drop the viscous term in the Burgers 
equation without giving some meaning to vPvP x at shocks. 
This can be done using BV- calculus ||, which allows one 
to write an equation for u° and gives rules for manipulat- 
ing the terms entering this equation and computing the 
effect of the viscous term in the inviscid limit. An alter- 
native, more intuitive, way of accessing the effect of the 
viscous shock on the velocity profile outside the shock 
is to carry out an asymptotic analysis near and inside 
the shock. Here we will take the second approach and 
refer the interested reader to |l(| for the first approach 
with BV-calculus. It is important to remark that the two 
approaches lead to the same results. 

Before considering velocity gradient, it is helpful to 
study the statistics of velocity itself. Let R(u; t) be the 
pdf of u[x, t). Assuming statistical homogeneity, R sat- 
isfies 



R t = B R UU - is((u xx \u}R) u , 



(7) 



where Bq = B(Q). To compute — vUu xx \u)R^ u , let us 
note that for v <C 1 , the solutions of (TO) consist of smooth 



1 



pieces where the viscous effect is negligible, separated by 
thin shock layers inside which the viscous effect is im- 
portant. Let u ut(x,t) be the solution of the Burgers 
equation outside the viscous shock layer; u ou t can be ob- 
tained as a series expansion in v. To leading order in v, 
M ou t satisfies Riemann's equation, ut + uu x = /. In order 
to deal with the shock layer, say at x = y, define 



V 



(8) 



and write v = vq + W\ + 0{v 2 ). To leading order, 
vo(z,t) satisfies (vo — u)vq z = vq zz , yielding vo(z,t) — 
u — (s/2) tanh(s2/4) where u — dy/dt and s is the jump 
across the shock. Consequently we have the following 
generic velocity profile inside the shock layer: 



u in (x,t) = u - 



s , / s(x — y) 

- tanh -i ?L 

2 \ \v 



0(1/). 



(9) 



The actual values of u and s are obtained from the match- 
ing conditions between u- ln and it ou t- In terms of v and 
the stretched variable z, they are 



lim v = lim u out =u±-. 



(10) 



It is well-known that s < 0. 

We will use (^|) to evaluate the viscous term in (R) . By 
definition pij, 



v(u xx \u)R — v lim — / dx u xx S[u — u(x,t)]. (11) 
L^-oo 2L J _l 

In the limit v — -> only small intervals around the shocks 
will contribute to the integral. So, we can split the inte- 
gral into small pieces involving only the shock layers and 
use the generic form of u- m in the layers to evaluate these 
integrals. To 0(v) : this gives 

v(u xx \u)R 



I™ ^fT?y2 f dx u inxx S[u- u in (x,t)} 

lL N j Jj-th layer 
/• r+oo 

p / dsduT(u,s;t) / dz vo zz S[u — vo(z,t)], 



(12) 



where in the second integral we picked any particu- 
lar shock layer and we went to the stretched variable 
z = (x — y)jv. Here N denotes the number of shocks 
in [—L,L], p = p(t) = limL^oo 7V/2L is the shock den- 
sity, and T{u, s; t) is the probability density of u(y, t) and 
s(y, t) conditional on the property that there is a shock 
at position y (T is independent of y because of statistical 
homogeneity). The last integral in (|lj) can of course be 
evaluated using the explicit form of vq. Another, more 
elegant, way to proceed is to use the equation for vq, 
(vq — u)vq z = vq zz , and change the integration variable 
from z to Vo using dzvo zz = dvoVo zz /vo z — dvo(vo — u). 
The result is 



lim v(u xx \u)R = ~p I ds 



«— s/2 



u+s/2 



du (u — u)T(u, s; t). 

(13) 



This equation gives an exact expression for the viscous 
contribution in the limit v — > in terms of certain sta- 
tistical quantities associated with the shocks. Of course, 
using (^||) in (0) does not lead to a closed equation since 
T remains to be specified. However, information can al- 
ready be obtained at this point without resorting to any 
closure assumption. For instance, using (|l3|) in (R) and 
taking the second moment of the resulting equation yields 
2B - 2e with 



1 



e= limi/(t£) = T^P<|*r>. 

v— >0 Iz 



(14) 



In particular, at steady state p(|s| 3 ) = 12B . 

Similar calculations can be carried out for multi-point 
pdf's and, in particular, for W(w; x,t), the pdf of the 
velocity difference w(x,z,t) = u(x + z,t) — u(z,t). It 
leads to an equation of the form 



W t = -WW X 



dw' W x (w';x,t) 



(15) 



+2[B - B{x)]W ww + H{w; x, t), 
where, to 0(x), H is given by 
H = p[wS{w;t) + (s)S{w)] 

+ 2p dw' S(w';t) -2p6(w) + 0(x). 



(16) 



Here 6(w) is the Heaviside function and S(s;t) = 
f duT(u,s;t) is the conditional pdf of s(y,t). By direct 
substitution it may be shown that the solution of ( |l5| ) is, 

to o(x 2 ), m 



W 



(1 - px)-Q (-; t) + pxS(w; t) + 0{x 2 ). (17) 
x V x J 



The first term in this expression contains Q(£; t), the pdf 
of the non-singular part of the velocity gradient, to be 



considered below (see (20)). This term accounts for those 
realizations of the flow where there is no shock in between 
z and x + z (an event of probability 1 — px + 0(x 2 )). 
This term also leads to the consistency constraint that 
lim-j^o W = S(w) (using lim^^o Q(w/x; t)/x = S(w)). 
The next term in (p"7j), pxS(w;t), accounts for the real- 
izations of the flow where there is a shock in between z 
and x + z (an event of probability px + 0(x 2 )). Equa- 
tion (^) can be used to compute the structure functions, 
(|w| a ) = J dw \ w\ a W. To leading order this gives 



(M») 



x a (\zn+o(x) 

xp{\s\ a ) + 0(x 



if < a < 1, 
if 1 < a, 



(18) 



where (|£| a ) = / d£|£| a Q- Using p(|s| 3 ) = 12B , we get 
Kolmogorov's relation for a — 3 



2 



12xB . 



(19) 



We now go back to the velocity gradient. Observe first 
that, in the limit v — > 0, the velocity gradient can be 
written as 



Notice that the exponential terms are irrelevant in these 
expression since s < 0. Equation (26) implies 



2u t s t 2f x 
hm v lz = T ± = ?±, 

z—>±oo S S S 



(27) 



,(x, t) = £(x, t) + s(yj)S(x — yj), (20) where the last equality is just the definition of £±. Note 

j that (27) can be rewritten as 

*t = -f (£-+£+), «t = !(£-- €+)+/*• (28) 

In the limit v — > these are the equations of motion 
along the shock. 

We can now evaluate the viscous contribution using 

v{Zx*\Z)Q = v Hm / £(*,*)]. (29) 

L—>oc ZL/ J 

The calculation is similar to the one for the velocity and 
eventually leads to 



where the yj's are the locations of the shocks, £ is the 
ncm- singular part of u K . Assuming homogeneity, a direct 
consequence of (pel) is 



(21) 



K) = (0 + p(s) = o. 



Unlike the viscous case where £ = u x , hence (£) =0, we 
have in the inviscid limit (£) = —p(s) ^ 0. Note also that 
the inviscid limit of the solutions of (^) converge to the 
pdf of £ only, which is still going to be denoted by Q. 

To evaluate F, there are two ways to proceed. One is to 
rewrite ( p"5| ) in terms of the pdf of (u(x + z, t) — u(z, t))/x 
and take the limit as x goes to zero. This is the approach 
taken in flQfl . The other is to evaluate (jj) directly. The 
two approaches amount to different orders of taking the 
limit x — » 0, v — » 0, and give the same result. Hence the 
two limiting processes commute. We will take the second 
approach and evaluate (Eh using the same basic idea as 
above. Here, however, we have to proceed more carefully 
with the shock layer analysis. Differentiation of (||) gives 



8u 



-sech 



2 ( s (x - y) 



0(1). 



(22) 



While the next order term in (Q) was negligible in the 
limit v — » 0, the O(l) contribution to £; n (a;,t) actually 
dominates the Olv^ 1 ) contribution at the border of the 
shock layer because the latter falls exponentially fast as 
the outer region is approached, whereas the former tends 
to constants, say, £±. In particular, the matching be- 
tween £ ut(a;,i) and £i a (x,t) involves the O(l) terms. To 
see how matching takes place, differentiating the expres- 
sion for u- m , we have £j n = v Vo z + v\ z + 0{v). The 
matching condition between £j n and £ ou t reads 

lim v\ z = lira £ aut = £±. (23) 

The equation for V\ is 

vat + (i>o ~ u)viz + viv 0z = v lzz + f X) (24) 

and, from the above argument, the only information we 
really need about vi is its values at the boundaries z — > 
±oo. Since vq z falls exponentially fast for large \z\, fl24| ) 
reduces to 

fit ± y ± 2 U iz = v izz + /*> 2->±oo, (25) 

where we used the asymptotic values of vq. Thus, as 

z — > ±oo, 



.2u t S t . 2f 



F&t) = ~ J dss[V.^,s;t)+V + ^,s;t)], 



(30) 



where V±(£, s; t) are the conditional pdf's of £±(y, i) and 
s(y, t). The appearance of £± in ( |30| ) is of course a direct 
result of the 0(1) term in @. 

We now use ( |30| ) in (|^) and analyze some consequences 

of 

Qt^^Q+{eQ) 5 + B 1 Q^ + F(C;t). (31) 
Taking the first moment of ([n]) leads to 

<0* = [£ 3 C» + f [<«£-> + <*£+>], (32) 

where we used f d£ £F = p[(s£-) + (s£+)]/2. On the 
other hand, averaging the first equation in (E8) gives 



(/»<-», 



[<*£-> + <*£+>] 



(33) 



«i ~ =p z 2 ± z + cf + cif e 



±p±S2/2 



(26) 



This equation uses the fact that shocks are created at 
zero amplitude, and shock strengths add up at collision. 
These are consequences of the fact that the forcing is 
smooth in space [ fL3| . Since (£) t = — (p(s))t from (pi]), 
the comparison between ( |32] ) and ( p3| ) tells us that the 
boundary term in ( [32] ) must be zero. Since Q > 0, £ 3 Q 
has different sign for large positive and large negative 
values of £. Therefore we must have lim£^ +00 £ 3 Q = 
and lim^-oo £ 3 Q = 0. This proves that Q goes to zero 
faster than |£|~ 3 as £ — > — oo and £ — > +oo. 

The analysis can be carried out one step further for 
the stationary case (Qt = 0). In this case, treating (|3l| ) 
as an inhomogeneous second order ordinary differential 
equation, we can write its general solution as Q — C\Q± + 
C2Q2+Q3, where C% and O2 are constants, Qi and Q2 are 
two linearly independent solutions of the homogeneous 
equation associated with ( |3l| ) , and Q 3 is some particular 
solution of this equation. One such particular solution is 



3 



* de me) 



Bi 



<e A 'G(C'), (34) 



where A = £ 3 /(3£i) and 



5 , ^(O 



(35) 



With this particular solution, it can be shown (see for 
details) that the realizability constraints imply that C\ — 
C*2 = 0, i.e. the only non-negative, integrable solution is 
Q = Q3. Furthermore, in order that Q actually be non- 
negative, F must satisfy 

> F > Ce 2 e"« 3/(3Sl) as ^ +00, (36) 

for some constant C < 0. Substituting into (0), we get 



C-\£\- 3 f^d? ?F(?) as^-c 



Q 



(37) 



as £ — ► +oc 
C~\£\- a with a > 3 as 



which confirms the result Q 
£ — > —00, and gives (3 = 1. 

The actual value of the exponent a depends on the 
asymptotic behavior of F. The latter can be obtained 
from further considerations on the dynamics of the shock 
( p8| ) . This is rather involved and will be left to . The 
result gives a = 7/2 which confirms the prediction of 
H . Here we will restrict ourselves to an interpretation of 
the current approach in terms of the geometric picture. 
Observe that the largest values of £± are achieved just 
after the shock formation. Assume that a shock is created 
at time t = 0, position x = 0, and with velocity u = 0. 
Then, locally 



x = ut — au + • • • . 

It follows that for t <C 1 the solutions of = ut 
u±, behave as 



u± 



a 



and solutions of 1 = £t — 3au 2 ^, behave as 

1 



2b 



(38) 
au 5 , 



(39) 



(40) 



Assuming that these give the dominant contribution to 
F(£) for large negative values of £, the asymptotic form 
of F is 

F~C dts(t){6[$-Z-(t)]+5[l;-!; + (t)}}, (41) 
Jo 

where C is some constant related to the statistics of the 



shock life-time and a, and s(i), £±(i) are given by (39) 
©. The evaluation of @ gives F - C|C|" 5/2 , and 
hence, 



Q ~ C_[^| 



-7/2 



Even though this argument gives only a lower bound for 
F at large negative values of £, further arguments pre- 
sented in Jl(| indicate that this lower bound is actually 
sharp. 
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